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LONG-TERM  GOALS 

The  objective  of  the  work  is  to  study  wave  evolution  in  cohesive  sedimentary  environments,  toward 
the  development  of  an  effective,  stochastic  model  for  wave  dissipation  in  these  environments.  The 
project  will  focus  on  prediction  of  surface  wave  evolution  over  relatively  large  spatial  scales  (scales  of 
100  wavelengths),  intermediate  depth  to  shallow  water. 

OBJECTIVES 

The  strong  dissipative  effects  cohesive  sedimentary  environments  have  on  waves  are  well  known,  but 
little  understood.  The  commonly  accepted  long  wave  paradigm  (only  low  frequency  motion  is  affected 
due  to  strong  interaction  with  the  bottom)  is  contradicted  by  observations  (Sheremet  and  Stone  2003) 
showing  strong  dissipation  also  in  high  frequency  bands,  where  the  direct  wave-bottom  interaction  is 
weak.  The  goal  of  this  project  is  to  investigate  short  wave  dissipation  and  develop  a  theoretical 
formulation  for  its  mechanisms,  and  is  the  results  to  develop  a  numerical  formulation  of  dissipation 
terms,  amenable  to  implementation  into  existing  stochastic  spectral  wave  models  (e.g.  SWAN). 

APPROACH 

The  observational  component  of  this  project  is  based  on  WAVCIS,  a  large  ocean  observing  systems  in 
operation  at  Louisiana  State  University.  The  system  is  operational  and  provides  comprehensive  real¬ 
time  in-situ  measurements.  The  project  enhances  WAVCIS’  capabilities  with  sediment  monitoring 
sensors. 

The  simulation  component  will  consist  of  a  fully  detailed  pilot  model,  integrated  in  a  quasi-operational, 
nowcast  mode,  with  the  observation  systems.  The  pilot  model  will  implement  theoretical  approaches 
derived  from  the  study,  and  will  serve  as  a  benchmark  for  skill  assessing  subsequent  parameterizations 
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of  wave  dissipation  developed  as  modules  for  operational  use.  The  resulting  simulation/observation 
system  (workbench)  will  help  identify  the  dominant  energy  exchange  mechanisms,  their  time/spatial 
scales,  and  ways  to  uncouple  the  surface  wave  and  fluid  mud/seabed  motions.  It  will  also  maximize  the 
exposure  of  development  ideas  to  field  data.  The  database  of  test  cases  will  be  established  during  the 
project. 

WORK  COMPLETED 

Recent  field  observations  show  unexpected  short-wave  sensitivity  to  sediment  type.  Sheremet  and 
Stone  (2003)  studied  the  effects  of  sediment  type  on  wave  propagation  by  comparing  synchronous 
wave  observations  from  two  locations  along  the  same  isobath  (5  m),  one  sandy  and  the  other  muddy. 
The  locations  were  on  a  very  mild  sloping  shelf  (slope  <  1/1000),  and  subject  to  nearly  identical 
meteorological  and  offshore  wave  conditions.  Comparable  short-wave  (f  >  0.2  Hz)  energy  levels  were 
recorded  at  the  two  sites  when  wind  forcing  was  significant.  However,  during  calmer  periods,  such  as 
in  the  wake  of  a  frontal  passage,  short  waves  dissipated  much  faster  at  the  muddy  site.  Understanding 
wave  dissipation  processes  in  muddy  environments  has  been  hampered  by  the  scarcity  of  field  data. 

Suspended  sediment  concentration  (SSC)  monitoring  devices  (OBS  -  Optical  Backscatter  Sensor)  have 
been  deployed  at  a  WAYCIS  station  located  on  a  muddy  shelf  (near  5-m  isobath).  To  understand  the 
formation  of  the  fluid-mud  layer,  SSC  data  were  analyzed  using  a  one-dimensional  numerical  model 
proposed  by  Li  and  Parchure  (1998)  to  simulate  a  non-equilibrium  suspension  profile.  The  focus  of 
simulations  was  to  describe  the  regime  dominated  by  settling  during  the  waning  stage  of  the  storm, 
when  the  gradual  decrease  in  wave-induced  turbulence  shuts  down  the  resuspension  processes.  The 
increase  in  concentration  in  the  fluid  mud  layer  suppresses  the  turbulence  further,  accelerating  the 
decay.  Net  bottom  resuspension  mass  flux  and  key  parameters  controlling  sediment  settling  were 
estimated  using  previously  published  data  (Kemp  and  Wells,  1987;  Wolanski,  1989;  Hwang,  1989). 

The  analysis  suggested  the  formation  of  a  sharp  lutocline  at  the  top  of  a  layer  with  concentrations 

higher  than  10  kg/m 3  and  thickness  about  1  m.  The  fluid-mud  layer  formation  coincided  with  strong 
dissipation  of  both  swell  and  high-frequency  (period  <  5s)  seas. 

The  short-wave  dissipation  mechanism  is  unknown.  Sheremet  et  al.  (2004)  hypothesize  that  nonlinear 
spectral  transfer  effectively  couples  short  waves  to  the  dissipative  muddy  bottom,  resulting  in  a 
spectrum- wide  influence  of  this  bottom.  To  investigate  this  hypothesis,  a  fully  detailed  unidirectional 
pilot  model  was  developed  which  incorporates  quadratic  nonlinear  wave  interactions  and  wave 
dissipation  effects  induced  by  the  presence  of  a  fluid  mud  layer.  Mud  induced  wave  dissipation  is 
modeled  using  the  approach  of  Jiang  and  Mehta  (1996).  The  wave  model  is  based  on  the  fully 
dispersive  deterministic  mild  slope  equation  developed  by  Agnon  et  al.  (1993),  and  Agnon  and 
Sheremet  (1997). 

A  major  difficulty  of  the  fully  detailed  model  is  the  necessity  to  solve  a  very  complicated  linear 
dispersion  relation  with  a  complex  wavenumbers  (e.g.  Dalrymple  et  al.  1991),  which  would  be 
inconvenient  in  the  context  of  a  numerical  wave  model.  As  a  first  attempt,  we  used  the  simplified 
model  of  Ng  (2000),  a  boundary  layer  asymptote  of  that  of  Dalrymple  and  Liu  (1978),  which  expresses 
the  wavenumber  as  a  perturbation  expansion,  allowing  complex  dissipative  wavenumbers  to  be 
functions  of  the  nondissipative,  real-valued  wavenumber.  This  mechanism  was  added  to  the  model  of 
Kaihatu  and  Kirby  (1995),  a  nonlinear  wave  model  for  simulating  the  triad  interactions  redolent  of 
waves  in  the  nearshore  and  surf  zones. 
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RESULTS 


Figure  1  shows  wave  dissipation  characteristics  based  on  typical  values  of  parameters  for  the 

Louisiana  shelf  (mud-layer  thickness  of  0.3  m,  water  depth  is  5  m,  and  mud  density  of  1,100  kg/m  ). 
Figure  2  shows  the  dissipative  effects  of  this  mud  layer  on  the  evolution  of  a  single  “difference”  wave 
triad  (frequencies  fj  =0.05  Hz,/j>=0.5  Hz,  and  fj=  0.55  Hz).  Long  wave  dissipation  is  very  strong  - 

amplitude  transmission  (fraction  of  initial  amplitude)  after  1,000  m  propagation  is  <1%.  The  short 
waves  are  barely  dissipated  (amplitude  transmission  factors  of  93%  and  99%  respectively).  Evolution 
with  nonlinear  3-wave  interactions  (solid  lines)  deviate  from  the  linear  evolution  (dashed  lines).  With 
high-frequency  modes  2  and  3  having  comparable  strong  transmission,  energy  is  pumped  from  the 
shortest  mode  (2,  magenta)  which  decays  faster  than  linearly,  into  the  middle  mode  (2,  blue),  which 
decays  slower  than  linearly  (Figure  2a).  Compared  with  linear  evolution,  the  low  frequency  mode  1  is 
sustained  for  a  longer  distance  (Figure  2b)  by  energy  transferred  from  the  shorter  modes,  consistent 
with  a  steady  energy  flux  directed  toward  long  waves.  Different  dissipation  rates  for  the  long  wave 
affect  the  short  waves  (not  shown),  while  the  nonlinear  rate  of  decay  of  the  long  wave  appears  to  be 
very  similar  -  slightly  slower  in  the  case  of  weaker  long  wave  dissipation. 

Numerical  results  for  a  full-spectrum  simulation  are  presented  in  Figure  3.  for  the  same  water/fluid- 
mud  conditions) .  The  spectrum  is  characteristic  of  locally  forced  short  wind  waves  in  deep  water,  a 
JONS  WAP  spectrum  with  a  peak  period  of  fpeak= 4  s,  and  significant  height  of  Hs^=  1  m.  Figure  4 

represents  the  evolution  of  the  energy  flux  spectrum,  normalized  by  the  initial  spectrum.  In  the  absence 
of  dissipation,  linear  evolution  conserves  modal  spectrum.  Figure  4a  shows  the  effect  of  dissipation  on 
the  evolution  of  the  spectral  flux.  Long  waves  are  strongly  attenuated  and  disappear  from  the  spectrum 
within  the  first  1  km  of  propagation.  Higher  frequencies  are  less  affected  by  the  presence  of  mud. 

Figure  4b  shows  the  nonlinear  (3-wave  interactions)  evolution  of  the  same  spectrum.  Mediated  by  the 
3 -wave  interaction  energy  exchange  mechanism,  the  strong  dissipation  in  the  long  wave  band  is 
draining  energy  from  modes  with / >  0.5  Hz  which  otherwise  would  experience  practically  no 
dissipation. 

In  parallel  to  the  development  of  a  fully  detailed  model  we  also  implemented  the  Ng  (2000)  model. 

The  simplified  boundary-layer  mechanism  was  added  to  the  model  of  Kaihatu  and  Kirby  (1995),  a 
nonlinear  wave  model  for  simulating  nonlinear  wave  propagation  in  the  shoaling  and  surf  zone.  Figure 
4  shows  numerical  simulations  of  the  evolution  of  a  narrow  banded  spectrum  (JONSWAP,  peak 
period  10  s,  significant  height  1  m)  of  waves  over  a  mildly  (1%)  sloping  beach.  The  figure  is 
comprised  of  two  false  color  plots  of  the  log  of  spectral  density,  with  water  depth  on  the  abscissa  and 
frequency  on  the  ordinate.  Figure  4a  shows  the  shoaling  of  the  spectrum  over  a  sandy  slope;  the 
amplification  of  the  second  harmonic  of  the  spectral  peak  are  apparent  as  the  shoreline  is  approached. 
Additionally  there  is  a  strong  spreading  of  the  energy  to  higher  frequencies  due  to  nonlinear  energy 
transfer.  Figure  4b  shows  the  evolution  of  the  same  spectrum  over  a  muddy  beach,  with  a  mud  layer  of 
1100  kg/m3.  In  contrast  to  the  sandy  case,  the  second  harmonic  amplification  occurs  much  earlier  in 
the  shoaling  process;  this  is  likely  due  to  boundary  effects  as  the  (detuned)  phase-locking  adjusts  from 
a  random  phase  distribution,  rather  than  nonlinear  triad  effects.  The  overall  energy  damping  is  greater 
in  the  mud  case,  and  the  equilibration  of  the  energy  transfer  to  fill  in  “valleys”  between  distinct  peaks 
in  the  spectrum  occurs  at  an  earlier  stage  in  shoaling. 
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IMPACT/APPLICATIONS 


Much  of  the  present  and  near-future  Navy  capability  on  predicting  regional  and  nearshore  processes 
assumes  a  sandy-non-cohesive  sedimentary  environment.  The  present  research  enhances  this  capability 
by  identifying  processes  and  developing  mechanisms  which  allow  expansion  into  areas  which  deviate 
from  this. 

RELATED  PROJECTS 

The  project  is  closely  related  to  the  ongoing  ONR  project  “Coupled  dynamics  of  waves  and  fluid  mud 
layers”  (PI  Sheremet,  UFL  and  Dr.  M.  Allison,  Tulane),  and  also  closely  aligned  with  the  large  multi¬ 
year  NRL  ARI  entitled  “Coastal  Dynamics  of  Heterogeneous  Sedimentary  Environments”  (PI.  Dr.  K. 
Todd  Holland). 
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Figure  1:  Linear  characteristics  of  wave  dissipation  induced  by  the  presence  of  a  fluid-mud  layer  (in 
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this  example  water  depth  is  5  m,  mud  layer  thickness  is  30  cm ,  mud  density  is  1100  kg/m  ).  a) 
wavelength  change  with  respect  to  a  rigid  (e.g.  sandy  bottom),  b)  amplitude  transmission  at  1  km 
(ratio  of  amplitude  after  1  km  propagation  to  initial  amplitude)  as  a  function  of frequency. 


1.5 


x 

M— 

1  s> 
0 
c 
0 

-O 

0 

N 

0.5  78 
E 

o 

c 


o 

0  10  20  30  40  0  10  20  30  40 


Figure  2:  Nonlinear  evolution  of  a  triad  over  a  fluid  mud  layer.  Water  depth  and  mud  layer 
thickness  and  density  are  the  same  as  in  Figure  1.  Frequencies  are  fj  =0.05  Hz  (red),  f 2=0. 5  Hz 

(blue),  and  f $=0.55  Hz  (magenta).  Transmission  factors  atx=l  km  are  (a,b)  (<  1%,  93%,  99%),  and 
(c,d)  (10%>,  93%,  99%)  for  fj,  and  ff  respectively.  Linear  long  wave  dissipation  is  more  than  ten 
times  larger  in  (a,b)  than  (c,d).  Dashed  lines  represent  the  linear  (i.e.  non-coupled  modes)  evolution. 
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Normalized  energy  flux 


Figure  3:  Simulated  evolution  of  a  short  wave  spectrum  JONSWAP  (peak  T=4  s,  significant  height 
H=1  m)  for  a  spectrum  over  a  fluid-mud  layer  (same  parameters  as  in  Figure  1).  Energy  flux 
density  is  normalized  by  initial  value,  a)  Linear  model  ( including  dissipation),  b)  Nonlinear  model. 
A  normalized  plot  of  amplitude  transmission  at  1  km  (Figure  lb)  is  also  shown. 


Figure  4:  Numerical  simulation  (Kaihatu  &  Kirby ,  1995;  Ng,  2000)  of  spectral  density  evolution 
over  a  mildly  sloping  beach  (log  scale  spectral  density  vs.  depth  and  frequency).  The  shoreline  is  on 
to  the  right  of  the  plots,  a)  A  sandy  beach,  a)  A  muddy  beach  (fluid  mud  density  1100  kg/m3). 
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